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1.  Introduction 

The  II-IV-V2  semiconductor  CdSiP2  with  uniaxial  crystal  structure  is  attractive  for  nonlinear 
optical  frequency  mixing  applications  because  of  its  high  nonlinearity  (d36  >  80  pm/V), 
adequately  large  birefringence  (-0.05),  along  with  wide  bandgap  energy  (>  2  eV)  [1,2].  Since 
the  thermal  conductivity  of  CdSiP2  (13.6  W/m/K)  is  relatively  modest  compared  to,  say  that 
of  GaP  (110  W/m/K),  it  is  important  to  lower  the  absorption  losses  in  the  material  to  obtain 
high  power  frequency  conversion.  One  way  to  obtain  lower  absorption  is  to  lower  the 
temperature  of  the  material.  Figure  1  shows  the  infrared  transmission  in  CdSilT  measured  at 
90  K,  295  K,  and  450  K  with  the  higher  transmission  at  low  temperature  indicating  possibly 
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lower  losses.  However,  accurate  low  temperature-dependent  refractive  index  values  for 
CdSiP2  are  not  readily  available  -  although  they  have  been  measured  at  room  temperature 
[3,4].  A  thermo-optic  dispersion  formula  for  CdSiP2  has  been  derived  [5]  over  a  temperature 
range  of  10  to  70  °C  (283  to  343  K)  and  a  wavelength  range  of  0.5  to  6.6  pm.  In  this  work  we 
present  refractive  indices  of  CdSiPi  measured  over  a  temperature  range  of  90  to  450  K  in  the 
wavelength  range  of  0.97  to  9.5  pm. 

2.  Measurement  technique 

The  method  used  for  refractive  index  measurement  was  outlined  by  Moss  [6]  and  used  more 
recently  for  accurate  measurement  of  refractive  indices  of  GaAs  [7].  A  single  crystal  of 
CdSiPn  grown  by  the  Bridgman  technique  was  polished  on  both  sides  and  thinned  to  160  ±  2 
pm.  The  crystal  was  XZ  oriented  -  that  is,  the  a-  and  c-axes,  respectively,  were  in  the  plane  of 
the  surface,  and  the  Y-axis  was  an  a-axis  perpendicular  to  the  surface.  Transmission  spectra 
of  the  samples  were  taken  using  a  PerkinElmer  FTIR  spectrometer  run  with  a  step  size  of 
0.125  cm-1  over  the  wavenumber  range  of  16,667  to  400  cm-1  (i.e.,  wavelength  range  of  0.6 
to  25  pm)  and  over  a  temperature  range  of  78  K  to  400  K.  Light  incident  on  the  sample  was 
polarized  using  two  wire  grid  polarizers  (Thorlabs  WP25M-UB  and  WP25H-K),  covering  the 
spectral  ranges  from  0.25  to  4  pm  and  3  to  30  pm,  respectively.  With  light  polarized  along  the 
crystal  X  (or  a)  axis,  the  ordinary  index  na  was  obtained,  and  when  light  was  polarized  along 
the  crystal  Z  (or  c)  axis,  the  extraordinary  index  ne  was  obtained. 


Fig.  1.  Temperature  dependent  FTIR  transmission  spectra  of  CdSiP2  for  light  polarized  along 
the  ordinary  axis. 

The  spectra  shown  in  Fig.  1  exhibited  a  series  of  fringes  having  extrema  occurring  at 
specific  wavelengths  at  which  the  light  transmitted  through  the  samples  interfered 
constructively  or  destructively,  if  the  sample  thickness  is  denoted  by  d,  then  constructive 
interference  at  a  wavelength  1  indicates  a  relationship: 

2n(A)d  =  mA  (1) 
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where  m  is  an  integer  at  the  fringe  maxima. 

From  the  sample  thickness  (d)  and  the  fringe  number  (in)  the  refractive  index  n  can  be 
determined.  Since  neither  parameter  can  be  directly  measured  with  sufficient  accuracy,  the 
refractive  index  values  determined  from  prism  measurements  [3]  were  used  to  determine  both 
the  values  d  and  the  room  temperature  fringe  number  at  wavelengths  near  1  pm.  From  the 
thickness  and  fringe  number,  the  index  values  at  longer  wavelengths  and  at  other 
temperatures  were  determined. 

3.  Results  and  discussion 


From  the  FTIR  spectra  of  the  160  pm  thick  wafer,  the  results  of  the  refractive  index 
measurements  over  a  temperature  range  of  90  K  to  450  K  and  wavelength  range  of  0.97  to  9.5 
pm  from  are  shown  in  Fig.  2  and  referred  to  as  the  ‘experimental  values’.  They  are  fitted  with 
a  Sellmeier  equation  of  the  form 


n 


2 


A(T)  + 


B(T ) 
A2-C 


+ 


D 

A2-E 


(2) 


and  the  fit  parameters  are  summarized  in  Table  1.  This  form  of  the  Sellmeier  equation  was 
chosen  to  be  consistent  with  the  (equivalent)  form  presented  in  Ref.  [5],  and  is  discussed  in 
Ref.  [8]  (p.  99,  Eq.  (38)). 

The  difference  between  the  ‘experimental  values’  shown  in  Fig.  2  and  the  values  obtained 
from  Eq.  (2)  is  the  fit  error,  which  is  plotted  in  Fig.  3.  There  are  possibly  various  origins  of 
this  error,  including  a)  the  presence  of  additional,  but  weak,  resonances  not  accounted  for,  b) 
the  particular  manner  in  which  temperature-dependence  was  included  (as  a  polynomial  in  the 
coefficients,  and  not  in  the  resonant  frequencies),  and  c)  minor  errors  in  finding  fringe  peaks 
from  the  spectra.  Nevertheless,  we  find  good  agreement  between  the  experimentally 
determined  phase  matching  angles  and  the  temperature  dependence  of  idler  wavelengths,  as 
will  be  shown  in  Section  4. 

To  our  knowledge  these  are  the  first  cryogenic  and  above  room  temperature  index 
measurements  for  CdSiP2. 


Fig.  2.  The  experimentally  obtained  values  for  the  ordinary  and  extraordinary  refractive  indices 
of  CdSiP2  at  different  temperatures. 
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2468  10  2468 

Wavelength  (pm) 


Fig.  3.  Fit  errors  between  Sellmeier  expressions  and  experimental  refractive  index  at  different 
temperatures,  for  both  the  ordinary  (left  column)  and  extraordinary  (right  column)  axes. 

Table  1.  Temperature-dependent  Sellmeier  fit  coefficients3 


Coef.  nQ  ne 

A  11.95  + 5.3479  x  llT1  T+  5.5894  x  10 -7  T2  1 1.438  +  5.5408  x  KT*  T+  5.0458  x  10 ~7  2* 

B  0.6134  +  9.4768  x  1(T5  T+  2.0148  x  l(T7 1  0.61584  +  3.8668  x  1(T5  T+  2.9901  x  1(T7  T1 

C  0.101733  0.11182 

D  2334.22  2021.26 

E_ 833.205_ 777.162 

aThe  temperature  T  is  in  Kelvin  and  wavelength  X  is  in  microns.  The  fits  are  valid  for  a  wavelength  range  of  0.97  to 
9.5  pm,  for  temperatures  between  90  and  450  K. 

4.  Comparison  of  predicted  and  experimental  parametric  conversion 
wavelengths 

Temperature  dependence  of  the  signal  and  idler  wavelengths  obtained  by  parametric 
frequency  down  conversion  of  pump  lasers  at  wavelengths  of  1.064  and  1.99  pm  has  been 
presented  in  the  literature  [9,  10].  In  Ref.  [9],  the  results  of  optical  parametric  generation 
(OPG)  of  a  (pump)  laser  at  1.064  pm  are  presented  for  a  5  mm  thick  CdSiP:  crystal  cut  at  a 
phase  matching  angle  6  =  90°  between  the  crystal  c  axis  and  the  propagation  direction  of  light 
through  the  crystal.  The  temperature  of  the  crystal  was  varied  from  -295  K  to  425  K  and  the 
wavelength  of  the  generated  signal  beam  was  measured. 

Using  the  temperature  dependent  Sellmeier  Equations  given  in  Eq.  (2),  the  signal  and  idler 
wavelengths  satisfying  the  phase  matching  condition  (for  6  =  90°,  with  the  pump  beam 
traveling  as  an  extraordinary  wave  and  the  two  generated  beams  traveling  as  ordinary  waves), 
were  determined  over  a  temperature  range  of  90  to  450  K.  These  results  are  shown  in  Fig.  4, 
along  with  the  experimental  results  from  Ref.  [9].  The  calculated  wavelengths  agree  very  well 
with  those  experimentally  measured. 
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Fig.  4.  Comparison  between  the  experimental  temperature-dependent  idler  and  signal 
wavelengths  [8]  (symbols)  and  the  predictions  using  temperature-dependent  Sellmeier 
equation  (solid  lines)  for  Type  I  phase  matching  at  critical  phase  matching  angle,  and  CdSiP2  is 
pumped  at  1.064  pm. 

In  the  first  optical  parametric  oscillator  (OPO)  demonstrated  in  CdSiP2  with  a  pump  laser 
at  1.99  pm  [10],  the  generated  wavelengths  were  measured  as  a  function  of  the  crystal 
orientation  (i.e.,  phase  matching  angle)  at  room  temperature,  or  about  295  K  (Fig.  5(a)).  They 
were  also  found  as  a  function  of  crystal  temperature  for  a  fixed  phase  matching  angle  (Fig. 
5(b)).  It  is  worth  noting  that  by  interpolating  the  experimental  data,  at  295  K  (the  temperature 
for  Fig.  5(a)),  the  idler  generated  in  Fig.  5(b)  is  4.38  pm;  from  the  experimental  data  in  Fig. 
5(a),  this  wavelength  is  generated  at  a  phase  matching  angle  of  46.84°  (a  little  different  from 
the  crystal  cut  of  0  =  46°  ).  Using  the  Sellmeier  expressions  reported  here  (Eq.  (2),  the  best  fit 
to  the  temperature -dependent  data  is  found  at  a  phase  matching  angle  of  47.04°.  If  the  0.2° 
difference  between  the  fit  angle  and  the  angle  found  from  the  data  is  attributable  to  a  small 
error  in  the  reported  phase  matching  angles,  then  the  Sellmeier  expressions  above  well  match 
not  only  the  temperature-dependent  data  but  also  the  angle-dependent  data.  (See  Fig.  5(a), 
where  the  calculated  curve  was  shifted  to  lower  angles  by  0.2°  to  correspond  to  the 
experimental  values  as  reported.)  The  relatively  small  discrepancy  between  the  theory  and 
experimental  results  can  be  attributed  to  experimental  errors  in  the  measurement  of  the  phase 
matching  angle  (as  discussed),  the  crystal  temperature  (either  in  Ref.  [10]  or  in  the 
measurements  reported  here),  or  in  our  assumption  of  295  K  for  room  temperature  (as 
recorded  in  Ref.  [10]). 
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Fig.  5.  (a).  Comparison  of  angular  dependent  between  the  experimental  results  (Ref  [9], 
symbols)  and  the  prediction  using  the  Sellmeier  equation  at  295  K  (solid  line),  shifted  by  -0.2° 

(see  text);  Fig.  5(b),  Comparison  of  temperature  dependent  signal  and  idler  wavelengths 
between  experimental  results  (symbols)  from  Ref  [9].  and  the  prediction  using  the  Sellmeier 
equation  (solid  lines). 

5.  Temperature  dependent  phase  matching  curves  at  different  pump 
wavelengths 

The  predicted  temperature-  and  angle-  dependent  phase  matching  curves  provide  guidance  as 
to  the  ideal  crystal  cut,  target  operating  temperature,  polarization  type,  pump  laser  wavelength 
selection,  etc.  All  of  these  must  be  a  part  of  the  design  of  a  %{~]  frequency  conversion 
experiment.  In  what  follows,  we  present  the  calculations  for  two  special  cases  of  interest,  as  a 
function  of  the  pump  wavelength.  The  more  general  phase  matching  curves  are  then  given  in 
Figs.  9-14  for  a  series  of  temperatures  and  laser  pump  wavelengths  in  the  region  in  which  the 
Sellmeier  fits  are  valid.  (In  these  latter  cases,  the  plots  show  phase  matching  out  to  10  pm, 
which  is  only  a  small  extrapolation  from  the  fit  limit  of  9.5  pm.) 

Non-critical  phase  matching  (i.e.,  0  =  90°)  is  desirable  in  most  cases  of  frequency  mixing 
because  the  generated  wavelengths  are  not  in  this  case  sensitively  dependent  on  the  crystal 
angle.  Non-critically  phased  matched  signal  and  idler  wavelengths  are  shown  in  Fig.  6  for 
three  temperatures  spanning  the  measured  range,  employing  each  of  the  three  phase  matching 
types  ( e  — >  oo  ,  e  — >  eo  ,  e  — »  oe )  which  result  in  phase  matching  curves  falling  within  the 
Sellmeier  equation's  region  of  validity.  Note  that  for  the  Type  II  phase  matching  curves 
(dashed  and  dotted  lines),  the  nonlinear  coefficient  is  zero  for  a  non-critically  phase  matched 
crystal.  In  order  to  use  those  curves,  one  must  Pine  away  from  0  =  90°;  this  will  recover  some 
finite  nonlinear  coefficient,  but  the  resulting  process  will  be  far  less  efficient  than  Type  I 
phase  matching. 
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Phase  Matching 

- e—>o  o 

---  e—>e  o 
.  e—>o  e 


Temperature 

-  90  K 

-  300  K 

-  450  K 


Fig.  6.  Non-critical  phase  matching  signal  and  idler  wavelengths  as  a  function  of  wavelength, 
temperature,  and  phase  matching  type. 


The  case  of  degenerate  optical  parametric  generation  (the  reverse  of  second  harmonic 
generation)  is  also  frequently  of  special  interest.  Figure  7  below  shows  the  angle  at  which  a 
CdSiPn  crystal  should  be  cut  in  order  to  produce  degenerate  OPG  when  pumped  at  a 
wavelength  between  1.4  and  4.8  pm,  at  a  temperature  between  90  K  and  450  K. 


Degenerate  Parametric  Down  Conversion 


Phase  Matching 

- e—>  o  o 

- e—>  o  e 

Temperature 

-  90  K 

-  300  K 

-  450  K 


Fig.  7.  Phase  matching  angle  for  degenerate  optical  parametric  generation  at  90  K,  300  K,  and 
450  K,  as  a  function  of  the  pump  laser  wavelength.  These  are  shown  for  the  two  phase 
matching  types  ( eeo  and  eoe )  for  which  phase  matching  occurs. 


Phase  matching  predictions  in  CdSilb  show  a  particularly  interesting  situation  when 
pumped  with  an  extraordinary  wave  near  2.3  pm.  At  this  pump  wavelength,  the  phase 
matching  curves  are  almost  perfectly  vertical:  at  a  fixed  angle  and  temperahire,  there  is 
simultaneous  phase  matching  over  a  bandwidth  of  about  3  pm,  depending  on  the  pump 


wavelength.  For  a  1  cm  thick  crystal  pumped  at  2.34  pm,  the  value  of  the  sine2 


curve 


is  above  0.99  over  an  idler  range  of  2.2  pm,  and  the  FHWM  is  3.1  pm.  There  is  therefore  a 
potential  to  generate  a  relatively  large  bandwidth  in  this  fashion,  although  it  is  sensitive  to 
small  (0.05°)  changes  in  angle  (see  Fig.  8(b)).  It  is  worth  noting  that  this  feature  of 
simultaneous  phase  matching  across  a  large  bandwidth  is  present  regardless  of  the  operating 
temperature;  at  a  different  temperature,  one  needs  only  go  to  a  slightly  different  crystal  angle. 
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2.34  jjm  pump:  Type  I  (e—>oo)  sine  (AkL  /  2)  L  =  l°cm 


Fig.  8.  (a).  Phase  matching  curves  for  CdSiP2  pumped  at  2.34  pm,  for  a  variety  of 
temperatures.  The  curves  are  notable  for  simultaneously  phase  matching  a  broad  range  of 
wavelengths;  Fig.  8(b).  Plot  of  the  sine2  function  for  a  1  cm  thick  CdSiP2  crystal  held  at  300  K, 
for  several  crystal  angles.  This  curve  is  a  useful  metric  for  the  nonlinear  conversion  efficiency. 

Crystals  for  optical  parametric  generation  are  most  commonly  pumped  at  commercially 
available  laser  lines.  Accordingly,  signal  and  idler  wavelengths  are  plotted  for  several  of  these 
pump  wavelengths  in  Figs.  9-14  below  using  each  of  the  phase  matching  types  that  enable 
critical  phase  matching  in  the  short-wave  and  mid-wave  infrared  spectral  regions.  These  are 
shown  as  a  continuous  function  of  crystal  angle,  for  temperatures  between  100  K  and  450  K. 


1.064  |jm  pump:  Type  I  [e—>oo)  1.064  pm  pump:  Type  II  (e  — *oe) 


Fig.  9.  Phase  matching  curves  for  CdSiP2  pumped  at  1.064  pm,  for  a  variety  of  temperatures. 
The  phase  matching  types  differ  between  the  two  figures. 


1.25  pm  pump:  Type  I  (e—>oo) 


1.25  pm  pump:  Type  II  (e—>oe) 


Fig.  10.  Phase  matching  curves  for  CdSiP2  pumped  at  1.25  pm,  for  a  variety  of  temperatures. 
The  phase  matching  types  differ  between  the  two  figures. 


8 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


Signal  Idler  (gm)  Signal  Idler  (gm)  Signal  Idler  (gm) 


1.34  |jm  pump:  Type  I  (< e —>o  o ) 


1.34  pm  pump:  Type  II  (e—>oe) 


Fig.  11.  Phase  matching  curves  for  CdSiP2  pumped  at  1.34  pm,  for  a  variety  of  temperatures. 
The  phase  matching  types  differ  between  the  two  figures. 


Fig.  12.  Phase  matching  curves  for  CdSiP2  pumped  at  1.55  pm,  for  a  variety  of  temperatures. 
The  phase  matching  types  differ  between  the  two  figures. 


2  pm  pump:  Type  I  (e  —>o  o) 


2  pm  pump:  Type  II  (e—>oe) 


Fig.  13.  Phase  matching  curves  for  CdSiP2  pumped  at  2.00  pm,  for  a  variety  of  temperatures. 
The  phase  matching  types  differ  between  the  two  figures. 
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4.6  pm  pump:  Type  I  (e  —>o  o) 


Temperature  Stability  for  4.6  pm  Pump 


Fig.  14.  (a).  Phase  matching  curves  for  CdSiP2  pumped  at  4.6  pm,  for  a  variety  of 
temperatures.  The  phase  matching  angle  attains  a  maximum  between  cryogenic  and  elevated 
temperatures,  as  seen  in  Fig.  14(b),  for  the  case  of  degenerate  optical  parametric  down 
conversion. 


In  the  case  of  a  4.6  pm  pump,  phase  matching  was  only  predicted  to  occur  for  Type  I 
(e  —>oo)  phase  matching  when  limiting  the  wavelengths  to  those  over  which  Eq.  (2)  is  valid. 
In  Fig.  14(b),  the  phase  matching  angle  is  seen  to  attain  a  maximum  near  250K.  Operation 
near  the  maximum  is  more  temperature  stable. 

6.  Summary 

A  temperature  dependent  Sellmeier  equation  for  CdSilT  was  obtained  over  a  wide  range  of 
temperatures  and  wavelengths.  A  good  match  was  obtained  between  the  temperature 
dependent  signal  and  idler  wavelengths  predicted  by  the  equation  and  the  experimentally 
measured  values  for  two  pump  wavelengths.  The  dependence  of  the  generated  wavelengths 
on  phase-matching  angle  and  temperature  was  determined  for  various  pump  wavelengths. 
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